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CN . Abstract 

^ . We derive the expression for the jerk parameter in f{R) gravity. We use the Palatini 

' variational principle and the field equations in the Einstein conformal gauge. For 

Q\ ■ the particular case f{R) = R — the predicted value of the jerk parameter agrees 

with the SNLS SNIa and X-ray galaxy cluster distance data but does not with the 
SNIa gold sample data. 
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bj)! 1 Introduction 

> 

'jjj . A particular class of alternative theories of gravity that has recently attracted 

5^ i a lot of interest is that of the f{R) gravity models, in which the gravitational 

Lagrangian is a function of the curvature scalar R [1] It has been shown that 
current cosmic acceleration may originate from the addition of a term R^^ to 
the Einstein-Hilbert Lagrangian R [2]. 

As in general relativity, f{R) gravity theories obtain the field equations by 
varying the total action for both the field and matter. In this work we use the 
metric-afiine (Palatini) variational principle, according to which the metric 
and connection are considered as geometrically independent quantities, and 
the action is varied with respect to both of them [3]. The other one is the metric 
(Einstein-Hilbert) variational principle, according to which the action is varied 
with respect to the metric tensor g^^,, and the affine connection coefficients 
are the Christoffel symbols of g^j^^. Both approaches give the same result only 
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if we use the standard Einstein-Hilbert action [4]. The field equations in the 
Palatini formalism are second-order differential equations, while for metric 
theories they are fourth-order. Another remarkable property of the metric- 
affine approach is that the field equations in vacuum reduce to the standard 
Einstein equations of general relativity with a cosmological constant [4]. 



One can show that f{R) theories of gravitation are conformally equivalent 
to the Einstein theory of the gravitational field interacting with additional 
matter fields, if the action for matter does not depend on connection [3,5]. 
This can be done by means of a Legendre transformation, which in classical 
mechanics replaces the Lagrangian of a mechanical system with the Helmholtz 
Lagrangian. For f{R) gravity, the scalar degree of freedom due to the occur- 
rence of nonlinear second-order terms in the Lagrangian is transformed into 
an auxiliary scalar field [5]. The set of variables {gij,u, 4>) is commonly called 
the Jordan conformal gauge. In the Jordan gauge, the connection is metric 
incompatible unless f{R) — R. The compatibility can be restored by a certain 
conformal transformation of the metric: g^^ — > /i^j^ = f'(R)gfj_i,. The new set 
{h^u, 0) is called the Einstein conformal gauge, and we will regard the metric 
in this gauge as physical. 



f{R) gravity models have been compared with cosmological observations by 
several authors [6,7] and the problem of viabihty of these models is still open 
(see [8] and references therein). Current SNIa observations provide the data 
on the time evolution of the deceleration parameter q in the form of q = q{z), 
where z is the redshift [9] . The extraction of the information from these data 
depends, however, on assumed parametrization of q{z) [10]. For small values 
of z such a dependence can be linear, q{z) — q^ + q-iz [9], but its validity should 
fail at 2; ~ 1. A convenient method to describe models close to AC DM is based 
on the cosmic jerk parameter j, a dimensionless third derivative of the scale 
factor with respect to the cosmic time [11,12]. A deceleration-to-acceleration 
transition occurs for models with a positive value of jo and negative qo. Flat 
ACDM models have a constant jerk j — 1. 



In this work we derive the general expression for the jerk parameter in f{R) 
gravity. We use the field equations in the Palatini formalism and the Einstein 
conformal gauge [13]. We find the current value of this parameter for the case 
f{R) — R — [2,7] and compare it with recent cosmological data [10]. 
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2 Palatini Vciriation in f{R) gravity 



The action for f{R) gravity in the original (Jordan) gauge with the metric g^i, 
is given by [13] 



Here, \/—gf{R) is a Lagrangian density that depends on the curvature scalar 
R — Rn,^{T p^)g'^'^ , Sm is the action for matter represented symbolically by 
and independent of the connection, and n — Tildes indicate quantities 
calculated in the Jordan gauge. 

Variation of the action Sj with respect to g^j^i, yields the field equations 

nm^. - \f{R)~g,. = kT,,, (2) 



where the dynamical energy-momentum tensor of matter is generated by the 
Jordan metric tensor: 

6S^ = ^J'd'x^T^jr: (3) 



and the prime denotes the derivative of a function with respect to its variable. 
From variation of Sj with respect to the connection T^^^ it follows that this 
connection is given by the Christoffel symbols of the conformally transformed 
metric [5] 

g,. = nkyg,.. (4) 



The metric g^i, defines the Einstein gauge, in which the connection is metric 
compatible. 

The action (1) is dynamically equivalent to the following Helmholtz action 
[5,13]: 

Sh^-^ [ d'x^l^g[f{4>{p))+p{R- m)] + S^{g,,,^), (5) 



where p is a new scalar variable. The function (j){p) is determined by 

df{R) 



dR 



= P- (6) 



3 



Prom Eqs. (4) and (6) it follows that 



ci>-Rn<i>), (7) 



where R = R^j^i^iV p^)g^" is the Riemannian curvature scalar of the metric g^^. 

In the Einstein gauge, the action (5) becomes the standard Einstein-Hilbert 
action of general relativity with an additional scalar field: 

Se^-^ I d'x^g[R-p-^{p)+p-^f{<l>{p))] + Sm{p-^g^.,i;). (8) 



Choosing (f) (which is the curvature scalar in the Jordan gauge) as the scalar 
variable leads to 

Se = -7^ I d'x^[R - 2Vm + S^{[f\<t^r'g,,,^l^), (9) 



where V{(f)) is the effective potential 

v(^) = ^/'W-W). (10) 



Variation of the action (9) with respect to g/^i, yields the equations of the 
gravitational field in the Einstein gauge [13]: 

1 kT 

Rp^i^ - 2^9i.u = jj^ - V{(l))gi,^, (11) 



while variation with respect to (f) reproduces (7). Eqs. (7) and (11) give 

0/'(0)-2/(0) = /.T/'(0), (12) 



from which we obtain (p = 0(T). Substituting (p into the field equations (11) 
leads to a relation between the Ricci tensor and the energy-momentum tensor. 
Such a relation is in general nonlinear and depends on the form of the function 

my 
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3 The jerk parameter in f{R) gravity 



The jerk parameter in cosmology is defined as [11,12] 
d 



where a is the cosmic scale factor, H is the Hubble parameter, and the dot de- 
notes differentiation with respect to the cosmic time. This parameter appears 
in the fourth term of a Taylor expansion of the scale factor around Oq: 



^^l + Ho{t- to) - \qoHl{t - toY + \3^Hl{t - to)' 

CLO ZD 

+ O[{t-to)% (14) 
where the subscript denotes the present value. We can rewrite Eq. (13) as 

j=q + 2q^-j^, (15) 



where q is the deceleration parameter. For a fiat AC DM model j — I [10]. ^ 

Prom the gravitational field equations (11) applied to a fiat Robertson- Walker 
universe with dust we can derive the 0-dependence of the Hubble parame- 
ter [13] 



and the deceleration parameter [7] 



We also have the expression for the time dependence of 0: [13] 

• ^ V6c{4>f'-2f)^<j)f-3f 

2f2 + 0/'/"-6//" ■ ^ > 



^ This identity can be easily verified from Eq. (15) for special cases where the 
deceleration parameter is constant: q = 1/2 (matter dominated universe) and q = 
— 1 (de Sitter universe). 
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Combining Eqs. (16-18) and using q — (j)q'{(f)) leads to 

q _ i8f(0f -2/)(0r-0/r-//o 

H (0/'-3/)2(2/'2 + 0/'/"-6//") ■ 



(19) 



From Eq. (15) we finally obtain 



j(0) = [26' r + - 750V'V/" - 120//'^ + 18 fr 

+1890/V7" - 162/V"] X m' - 3f)\2f" + 4>ff" - 6//")]-'- (20) 

Wc now examine the case f{R) = R — f^, where a is a constant, which 
is a possible explanation of current cosmic acceleration [2]. In this model the 
present value of is 0o = (-1.05±0.01)a, where a = {7.35l:\:u)xl0--^^m~'^[7]. 
We do not need to know the exact value of a since it does not affect non- 
dimensional cosmological parameters. Substituting 0o into (20) gives 

Jo = 1.011°:°?. (21) 



This value does not overlap with the value j = 2.16^.^^^, obtained from the 
combination of three kinematical data sets: the gold sample of type la super- 
novae [9], the SNIa data from the SNLS project [14], and the X-ray galaxy 
cluster distance measurements [10]. The origin of this disagreement could come 
from the assumption of constant jerk used there. However, two of the three 
data sets separately are consistent with the f{R) = R ~ model: the SNLS 
SNIa set gives j = 1.32l};2i and the cluster set gives j = O.Sll^oo' it is 
the gold sample data that yields larger j = 2.75t['^l [10]- ^ 

In the f{R) — R — ^ model the deceleration-to-acceleration transition oc- 
curred at 0t = —^JbJSa [7]. The cosmic jerk parameter at this moment can 
be found from Eq. (20): 



This value shows that the jerk parameter in f{R) gravity changes significantly 
between the deceleration-to-acceleration transition and now, indicating the de- 
parture of f{R) gravity models from AC DM. It would be interesting to gen- 
eralize the kinematical approach of [10] to time dependent jerk and compare 

^ The value go = —0.81 it 0.14 found in [10] from the combined three data sets 
agrees with qq = — 0.671qq3 derived in the f{R) = R — ^ model [7]. Each set 
separately agrees with our model as well. 
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the results with f{R) gravity models. More constraints on these models could 
also be provided by non-dimensional parameters containing higher derivatives 
of the scale factor, such as the snap parameter s = [12]. 



4 Summary 



We derived the expression for the cosmic jerk parameter in f{R) gravity for- 
mulated in the Einstein gauge. We used the Palatini variational principle to 
obtain the field equations and apply them to a flat, homogeneous, and isotropic 
universe filled with dust. The value of the jerk parameter for the particular 
case f{R) = R — does not overlap with the value obtained from cosmo log- 
ical data of the SNIa gold sample, but is consistent with the values obtained 
from more recent SNLS SNIa data and the X-ray galaxy cluster data [10]. 
Therefore, Palatini f{R) models in the Einstein gauge, including the case 

2 

f{R) = R — f^, provide a possible explanation of current cosmic accelera- 
tion. Further observations should give stronger constraints on j and on f{R) 
gravity. 
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